Realization of Discrete Quantum Billiard in 2D Optical Lattices 



Dmitry O. Krimer^'^ and Ramaz Khomeriki^'^ 

^ Theoretische Physik, Universitdt Tubingen, Auf der Morgenstelle I4, 72076 Tubingen, Germany 
^Max-Planck Institute for the Physics of Complex Systems, Nothnitzer Str. 38, 01187 Dresden, Germany 
^Physics Department, Tbilisi State University, Ghavchavadze 3, 0128 Tbilisi, Georgia 

We propose the method for optical visualization of Bose-Hubbard model with two interacting 
bosons in the form of two-dimensional (2D) optical lattices consisting of optical waveguides, where 
the waveguides at the diagonal are characterized by different refractive index than others elsewhere, 
modeling the boson-boson interaction. We study the light intensity distribution function averaged 
over direction of propagation for both ordered and disordered cases, exploring sensitivity of the 
averaged picture with respect to the beam injection position. For our finite systems the resulting 
patterns reminiscent the ones set in billiards and therefore we introduce a definition of discrete 
quantum billiard discussing the possible relevance to its well established continuous counterpart. 
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A very rich variety of wave phenomena that have orig- 
inally been discovered in the context of atomic and solid 
state physics attracted recently much attention due to 
their deep analogy to optical systems. A first prominent 
example is the Anderson localization, the phenomenon 
which was originally discovered as the localization of elec- 
tronic wavefunction in disordered crystals [l^ and later 
understood as a fundamental concept being universal 
phenomenon of wave physics. Related recent experiments 
were performed on light propagation in spatially random 
nonlinear optical media [2, 3] and on Bose-Einstein con- 
densate expansions in random optical potentials Q. A 
second example is the well known solid state problem 
of an electron in a periodic potential with an additional 
electric field, which lead to investigations of Bloch os- 
cillations and Landau- Zener tunneling in various phys- 
ical systems such as ultracold atoms in optical lattices 
[il-Q and optical waves in photonic lattices [s', Q . A re- 
cent progress in the experiments stimulated a new turn 
in theoretical studies dealing with the evolution of a wave 
packet in nonlinear disordered chains [10], in a nonlinear 
Stark ladder [11] and the effect of Anderson localization 
of light near boundaries of disordered photonic lattices 
p!2l| which are just a few recent examples. A third very 
interesting example is a classical analog of beam dynam- 
ics in one-dimensional (ID) photonic lattices to quantum 
coherent and displaced Fock states [13] and a classical 
realization of the two-site Bose-Hubbard model (appli- 
cable to the physics of strongly interacting many-body 
systems), based on light transport in engineered optical 
waveguide lattices [l^ . 

In this Letter we study a classical analog of beam prop- 
agating in 2D photonic lattices to quantum coherent dy- 
namics of two particles in one-dimensional chain using 
the Bose-Hubbard model. We consider different situa- 
tions ranging from the simple ordered case without inter- 



action to the disordered case with interaction in our finite 
systems. Sometimes the resulting patterns look pretty 
similar as the ones for the classical and/or quantum bil- 
liards which are known to exhibit regular and chaotic 
behaviors (see e.g. Ref. f\3\). We would like to empha- 
size particularly the growing interest to the two-particle 
problem in the context of quantum correlations between 
two noninteracting particles evolving simultaneously in 
a disordered medium [l6| and quantum walks of corre- 
lated photons which provides a route to universal quan- 
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FIG. 1: Geometry of setup: A beam enters into 2D optical lat- 
tice, and propagates along the z axis. The refract ivity index 
is invariant along the z axis and either periodic or disordered 
in transverse directions. The corresponding mapping to the 
dynamics of two interacting distinguishable bosons in a chain 
is also done (see the text for details) . The interaction between 
bosons is introduced by taking the refractive index for the di- 
agonal waveguides different from the rest. The injection of 
a beam to the diagonal waveguide mimics lunching initially 
both bosons at the same site (upper inset), while injecting 
the beam into off-diagonal waveguide corresponds to the two 
bosons located initially on different sites (lower inset). 



FIG. 2: (Color online) Characteristic pictures for discrete 
quantum billiard realization for different injection points de- 
picted by a red cross in the absence of disorder and interaction 
(W = and [7 = 0). In the main graphs the averaged PDFs 
[see Eq. ([3])] are displayed and accompanied by the lower in- 
sets which show the numbers of initially injected waveguides 
(or number of sites at which the particles are initially located) . 

turn computation [l^. Thus, the obtained results might 
be apphcable to both classical and quantum systems. 

Let us introduce a standard Bose-Hubbard Hamil- 
tonian describing two distinguishable bosons (or two 
fermions with opposite spins) in a chain with N sites 

N 

where 6^ (a^) and hj (aj) are bosons creation and an- 
nihilation operators on a lattice site j and U is the on- 
site interaction strength. Starting from the time depen- 
dent Schrodinger equation idt\'^{t)) = 1-L\^{t)) we ex- 
pand \^(t)) in terms of the N'^ orthonormal eigenstates 
of a number operator, |m,n) = 6+a+|0), as \^(t)) — 

N 

Cmn{t)\m^n), where the amplitudes Cmn{t) satisfy 

m,n=l 

to the following set of equations 

N 

m',n' = l 

Note that Eq. ([2j) is invariant under permutation of m 
and n, and therefore it is natural to represent Cmn as 
a sum of symmetric = {cmn + Cnm)/V^ and anti- 
symmetric = {Cmn — Cnm)/V^ fuUCtioUS. In SUCh 



FIG. 3: (Color online) All parameters and quantities are the 
same as in Fig. [21 except the interaction constant t/ = 1. 



a basis the matrix R^r^ decomposed into two irre- 
ducible parts, one of which corresponds to the Bose- 
Hubbard model with two indistinguishable bosons and 
the other describes the physics of two indistinguishable 
spinless fermions. For the symmetric initial conditions, 
Cmn(O) = Cnm(0)7 the dynamics is reduced to the former 
case (two indistinguishable bosons on sites m and n), 
whereas the latter case is realized for the antisymmetric 
initial conditions, Cmn(O) = — Cnm(O). 

Remarkably, Eq. (|2]) is the same as one used for the 
description of light propagation through 2D optical lat- 
tices (see Fig. [T]) within the tight-binding approximation, 
where longitudinal dimension z plays a role of time. This 
approximation is valid when a lattice is constructed such 
that tunneling into nearest neighboring waveguides is al- 
lowed only and there is a difference between the refrac- 
tive indices of the diagonal Ud and off-diagonal no waveg- 
uides which models the interaction (with the interaction 
strength U ^ uq — rid). Thus, injecting light beam at 
the waveguide with a position x = y = n (asym- 
metric initial conditions) corresponds to the dynamics of 
two distinguisable interacting bosons in a chain, placed 
initially on sites m and n. One can also think about 
the Bose-Einstein condensate embedded into 2D optical 
lattice and then Eq. (|2]) describes the evolution of some 
initial matter wave packet through the lattice. 

In this Letter we consider the system with hard bound- 
aries having Cmn = outside a square and monitor the 
time averaged wave function 

1 

Pmn = lim — / \Cmn{t)fdt, (3) 



3 




FIG. 4: (Color online) Characteristic pictures for discrete 
quantum billiard realization with a rigid square placed inside 
a system for the interaction constant U — 1 and the same in- 
jection point depicted by a red cross. (a),(b): symmetric and 
asymmetric situations, respectively, with the corresponding 
probability density functions of eigenvalue spacings s shown 
in (c) and (d). (c): the Poisson distribution JH]), with the av- 
erage spacing d — 0.0085. (d): the Wigner-Dyson distribution 
([7|), with the average spacing d — 0.0075. 



referring to Pmn as to the averaged two-particle probabil- 
ity distribution function (PDF). To calculate PDFs we, 
at first, solve the eigenvalue problem i-L\q) = Wq) and 
then expand Cmn{t) with respect to the eigenvectors as 



-iX„t 



(4) 



q=l 



where C 



{q\m^n) is the eigenvector which belongs 



to the eigenvalue Xq and 



- X/m n=l ^rnn{^)^'mn IS itS 

initial amplitude. Next, the averaged PDF is calculated 
by the following formula: 



(5) 



where the first sum runs over all nondegenerate eigen- 
values, whereas the second sum corresponds to the sum- 
mation with respect to r-fold degenerate eigenvalues Xq. . 

Intuitively it seems that the light injected into one of 
the waveguides should spread over a whole lattice, how- 
ever the real situation is completely opposite due to the 
interference from the hard boundaries. Let us start form 



the simplest noninter acting case, U = (for the optical 
counterpart shown in Fig. [1] waveguides must all be iden- 
tical). As is seen from Fig. [2] a well defined pattern for 
Pmn corresponds to each initial injection point. Thus, 
the system keeps the information about its initial state 
and from the averaged picture one can recover an initial 
signal. It should be noted, that these patterns might be 
strongly modified when the interaction is switched on, 
U ^ (see Fig. [3j). Remarkably, the pattern's structure 
reminiscent the one sets in billiards, therefore, we intro- 
duce the notion of discrete quantum billiard and seek for 
the analogies with usual continuous counterparts. The 
first step towards this direction is to explore the possi- 
bility of quantum chaos realization in such systems. We 
consider two possibilities to observe the transition to- 
wards quantum chaos. The first one is a symmetry break- 
ing by placing a square with rigid boundaries inside the 
system as is shown in Fig. HJ We monitor then the statis- 
tical properties of the eigenvalue spacings s = |Aq+i — Xq\ 
for different locations of the square, keeping injection 
point and interaction constant the same. It is seen that 
in the symmetric case the Poisson distribution 



P{s) = 1/d-e 



-s/d 



(6) 



is realized, while for the asymmetric case the Wigner- 
Dyson distribution is observed 



P{s) = 7Ts/{2d^) • e- 



(7) 



Thus the onset of quantum chaos could be visualized via 
the classical optical system of coupled waveguides. 

The second mechanism of quantum chaos realization is 
an introduction of the disorder via adding the following 
terms to the Bose-Hubbard Hamiltonian ([1]) 



N 



(8) 



where for the sake of simplicity we take symmetric dis- 
order = = Cj {cj are random numbers from the 
interval [—W/2^W/2]^ W being the disorder strength). 
In the optical context presented in Fig. [1] it implies the 
usage of symmetric (under permutation of m and n) dis- 
ordered lattice described by the following modified evo- 
lution equations: 



N 

E 



(9) 



where the matrix R^n 



is the same as in Eq. (|2]) and 
Wmn=^m + are Correlated disorder parameters. An 
uncorrelated disorder can also be taken with Wmn being 
a sum of two independent random numbers for each m 
and n. In the optical context such a situation might be re- 
alized by taking either correlated or uncorrelated random 
refractive index distribution in a whole lattice. Note, that 
only the former case coincides with the dynamics of two 
bosons in a disordered chain with a random potential. 
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FIG. 5: (Color online) Averaged PDFs both in time and over 
many disorder realizations (all parameters are shown on the 
figures). Injection point in all graphs is taken at the middle 
of 2D optical lattice. 




We consider two particles which are initially launched 
on the same site at the middle of a chain, m = n = N/2 
(the optical counterpart corresponds to the beam injec- 
tion at the central waveguide in 2D optical lattices). The 
typical structures for Pmn^ averaged out with respect to 
many disorder realizations, are shown in Fig. O As is 
seen, the averaged PDFs demonstrate well pronounced 
patterns, which look differently as compared to the case 
with a single disorder realization, when PDF has many 
spots at different locations. For the noninter acting case 
PDF has an anisotropic structure with two distinct direc- 
tions, m = N/2 and n = N/2, along which the particle 
motion mostly develops in average. For the case with 
correlated disorder, a slight contribution of two particle 
states is also visible. For [7 = 2, the interaction is already 
strong enough, such that the contribution of states cor- 
responding to the breather band becomes essential and 
the two particles mostly prefer to form a composite state 
and travel together. Interestingly, the level of system's 
chaoticity might easily be governed by setting different 
values of disorder and interaction strengths. To demon- 
strate this we study the statistical properties of the eigen- 
value spacings s = \Xq-\-i — Xq\ for three different cases as 
is shown in Fig. [6] (the events with 5 = due to degener- 
acy are not counted). It is seen that for nonzero disorder 
and interaction strengths, probability density function of 
the spacings s has a tendency to go to the Wigner-Dyson 
distribution and, as a consequence, a system becomes 
more chaotic. 
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Concluding, in this paper we have discussed various 
interpretation of the optical beam propagation through 
2D optical crystal ranging from interacting cold atom 
dynamics and two particle Anderson localization to the 
quantum billiard problems connected with transition to 
quantum chaoticity. 



FIG. 6: (Color online) Probability density functions of eigen- 
value spacings s for three different cases and a chain with 
iV = 51. Black curve: W=0, U=2 (disorder strength is zero 
only). Red curve (dark gray): W=2, U=2. The case with un- 
correlated disorder is considered. Green dash-dotted curve: 
the Wigner-Dyson distribution dT)), with the average spacing 
d — 0.0032. Blue dashed curve: the Poisson distribution (|6|), 
with the average spacing d — 0.0032. 
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